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O ■ Abstract 
O 



This paper shows that an integrable approximation of the spring pendulum, when tuned 
to be in 1 : 1 : 2 resonance, has monodromy. The stepwise precession angle of the swing 
■ plane of the resonant spring pendulum is shown to be a rotation number of the integrable 

CN| . approximation. Due to the monodromy, this rotation number is not a globally defined 

function of the integrals. In fact at lowest order it is given by arg(a + ib) where a and 
b are functions of the integrals. The resonant swing spring is therefore a system where 
' monodromy has easily observed physical consequences. 

ri : 

1 Introduction 

The spring pendulum or swing spring is one of the simplest possible mechani- 
qq \ cal systems. It is a spring with one end fixed and at the other end is attached 

a mass that is acted on by a constant vertical gravitation field. The name 
swing spring comes from the fact that for appropriate initial conditions the 
t— i . mass can either swing like a pendulum or bounce up and down like a spring. 

However, if the frequencies of the swinging and springing motion in linear 
approximation near the equilibrium are in resonance, these two types of mo- 
tions are intricately intertwined. In particular, the following motion is easily 
observed. Starting with the weakly unstable vertical springing motion, the 
. system evolves into a planar swinging motion. This swinging motion is also 

transient and the system returns to its original springing motion. This cycle 
then repeats. Lynch observed that the orientation of the swing plane 
typically changes from one swinging phase to the next. Moreover, the angle 
between the swing planes of any two successive swinging phases is constant. 
However, the angle between the swing planes depends on initial conditions. 
He called this phenomenon the stepwise precession of the swing plane of 
the swing spring. It is this phenomenon that we are going to explain both 
qualitatively and quantitatively. 
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The swing spring has a long history that is well described in J7[ . The ear- 



liest comprehensive work on the planar spring pendulum is [TT]] . This paper 
gives a classical treatment of the 1 : 2 resonance using action angle variables. 
It is written in the spirit of the old quantum mechanics and was actually 
motivated by the Fermi 1 : 2 resonance in CO2. The advent of modern 
quantum mechanics seems to have made this type of analysis old-fashioned 
if not archaic. As Lynch points out most of the previous work is only 
concerned with the planar spring pendulum, so that the stepwise precession 
of the swing plane cannot be found. Some progress on the three dimensional 
system was made in pi. After that Holm and Lynch || found that the sys- 
tem can be approximated by the 3-waves system and derived a differential 
equation for the angle of the swing plane. This was done using "pattern 
evocation in shape space" ||. We show that equation found by Holm and 
Lynch is exact and is nothing but the equation for the evolution of one the 
angles of the action-angle coordinates of an integrable approximation to the 
resonant swing spring. We then trace the origin of the stepwise precession 
to the existence of Hamiltonian monodromy in this integrable approxima- 
tion. Hamiltonian monodromy is an obstruction to the existence of global 
action variables, which was first described in || (see also 0). It generically 
appears around an equilibrium point of an integrable two degree of freedom 
Hamiltonian system whose linearization has a complex quartuple ±a ± % /3 



of eigenvalues fll2[ . Such an equilibrium point is of focus- focus type. The 
integrable approximation of the resonant swing spring has three degrees of 
freedom. But after reduction of a symmetry, one obtains a two degree of free- 
dom system with a focus-focus point as a relative equilibrium. Physically this 
corresponds to the pure springing motion of the system. For the purpose of 
the present paper the most important consequence of monodromy is that the 
rotation number of invariant tori, that is, the ratio of their frequencies, near 
the singularity is not a single valued function. Our main result is that the 
stepwise precession of the swing plane is given by such a rotation number, 
which explicitly has the form 



A$ = arg(a + ib) + 0(Va 2 + b 2 ). (1) 

Here a and b are simple function of the integrals of the system, and a, b — > 
near the equilibrium of the swing spring. The amazing feature of ([I]) is that 
it is not differentiable at the origin. This means that no matter how small 
the initial perturbation from the equilibrium is, one can always obtain all 
possible values for A$. In general, such multivaluednes has been described 



2 



for integrable foliations near focus- focus singularitites by Vu Ngoc [ID]. In 
some sense our result is a special case of his. However, he did not study the 
influence of the Hamiltonian, but only the foliation. 



The paper is organised as follows. In section 2 we briefly recall the physics 
of the swing spring. Then we derive an integrable approximation which is 
valid near the resonant equilibrium point. This integrable system is then 
reduced to a one degree of freedom system in section 4. We then describe the 
geometry of the image of its energy momentum map and show that there is 
monodromy. The dynamics of the swing angle $ is described in section 6 and 
finally we obtain equation (|1]) for the rotation number Afi, by approximating 
an elliptic integral. 



2 The physics of the swing spring 

The spring pendulum is a point particle r = (x, y, z) in M 3 of mass m attached 
to a spring which moves in a constant vertical gravitation field. Its potential 
energy is 

V(r)=mgz + lk(£ -\\r\\) 2 , (2) 

where ||r|| = ^/ x 2 + y 2 + z 2 . The unstretched length of the spring with 
spring constant k is £q and g is the constant of gravity. The motion of the 
swing spring is governed by Newton's equations 

mf = — grad V(r). (3) 

The system is in equilibrium when the forces of gravity and the spring bal- 
ance, that is, when 

grady(r) = 0. (4) 

Thus x = y = and z = —£. The equilibrium length £ of the spring is 
determined by 

k{£ - 4) = mg. (5) 

Measuring length in units of equilibrium length £, mass in units of m 
and time in units of (which is the period of small amplitude pendulum 
oscillations), we see that the Hamiltonian of the swing spring on phase space 
T*M 3 with canonical coordinates (x,y, z,p x ,p y ,p z ) is 

H =l{p 2 x +p 2 y +p 2 z ) + U{x,y,z), 
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where 




and v 




Since t > £q, we have v > 1. This says that the 



frequency of the spring oscillation is greater than the frequency of the small 
amplitude pendulum oscillations. The reason for this is that the frequen- 
cies are not independent because they are coupled by (|). Expanding the 
potential U @ about its stable equilibrium (0,0, —1) to cubic terms gives 



cubic approximation to the Hamiltonian of the swing spring in dimensionless 
variables with its equilibrium shifted to the origin is the Hamiltonian 



on phase space T*M. 3 with canonical coordinates Z = (x,y, z,p x ,p y ,p z ). By 
rescaling the coordinates and changing the time scale, we may consider p, to 
be a small parameter, which measures the distance to the origin. 

The swing spring Hamiltonian undergoes a 1 : 1 : 2 resonance when v = 2. 
In the remainder of this paper we will study only this case. The Hamiltonian 
of the cubic approximation of the resonant swing spring is 



H = H + V = l( P l+ P l+ P l) + l( x ^+y 2 + 4,z 2 )-fi{x 2 + y 2 )z. (9) 

3 An approximating integrable system 



In this section we find an integrable approximation to the resonant swing 
spring (TC, T*M 3 , uj = dx A dp x + dy A dp y + dz A dp z ) by averaging over the 
flow of the quadratic part of 7i. Due to the resonance there will be secular 
terms at cubic order. From general results it is known that a large measure of 
initial conditions of the averaged system stays close to the original solutions 
for long times, see e.g. [[![]. 

Clearly the Hamiltonian 7i (^) is invariant under rotation about the z-axis 




H u = \{pl+pl+ pi) + | (x 2 + y 1 + z/ 2 z 2 ) - M* 2 + V 2 )z 





















ip : (s,Z) h-> <p a (Z) 
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which has momentum 



xp y - yp x - 



(11) 



The quadratic terms H of the Hamiltonian 7i (0) are in 1 : 1 : 2 resonance. 
The flow of the Hamiltonian vector field Xh corresponding to H Q generates 
the oscillator action 
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(12) 



To make the oscillator action a symmetry, we average the Hamiltonian Ti, 
over the integral curves of Xh - We obtain 



= H n 



2t 



Hoip t {Z) dt 

| [{xPx + yp y )Pz + {x 2 + y 2 )z - (p 2 x + p 2 y )z] . (13) 



By construction the Hamiltonian 7i is invariant under the oscillator action 
(|T2l). Therefore {7i, Hq} = 0, where {, } is the standard Poisson bracket 
corresponding to to. It is straightforward to check that H is also invariant 
under the axial action (ITof) and hence {7i, L} = 0. Since the axial and 
oscillator actions commute, it follows that {Hq, L} = 0. Consequently, the 
system (TC, H , L,T*R 3 ,u) is Liouville integrable. 

For later purposes it is useful to study a different, but equivalent Liouville 
integrable system. This equivalent system is obtained by a linear symplectic 
transformation that diagonalizes both, H and L. Consider the invertible 
linear map 



* : T*R 3 -> T*IR 3 : (f , 77, C, P{, P„ P C ) ^ ( x > V> z > P*> Py> P*) = 

= ^(Pv + €iPt + VX,P£-V,Pv-£i 2 Pc) ■ 
Then ^ is symplectic, that is, 

\$j*(uj) = d£A dp;: + drj A dp v + d( A dp ? = w. 
Moreover, \I/ diagonalizes the momenta L (|TT| ) and i^o ®, namely, 
L = **(L) = W n + r, 2 )-\{p\ + e) 



(14) 



(15) 
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and 



H = V(H ) = l( P l + e+p 2 rj + V 2 + 2(pl + C 2 )) (16) 
In new coordinates H = (£, T], CiPtiPwPt) 011 (T^^ 3 ,^) the averaged Hamil- 



tonian H flIS|) becomes 

H = ^* (H) = H + \ [(Cp c - CPt)v - (£C + P#<)p,] , (17) 
where A = Hamilton's equations for the integral curves of Xg are 



£ =Pe-A« + Pi,Pc) 
77 =p„ - A(£C + p ? Pf) 
C = 2p ( + \(£r] - 
Introduce new momenta 

J 1 = § (#0 - L) 



p ? = -£ - A(r/p c - Cp^) 

Pr? =-v- H&c ~ Cp?) ( 18 ) 
Pc = -2( + H&r, + vPt)- 



wi+e+Pi+e 



(19) 



and 



J 2 = \{H + L) = l^ + ^+pl + C 2 ), 



(20) 

whose Hamiltonian vectors fields Xji and X.p have flows giving the S 1 - 
actions 



¥>f = (*,S) 



and 
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(21) 



(22) 



respectively. Since {J 1 , J 2 } 
leave the Hamiltonian H (|T~ 



= 0, these actions commute. Moreover, they 
invariant. Thus (H, J 1 , J 2 , T*R 3 



uj) is a Liou- 
ville integrable system. 

Equation (|18D is a variant of the 3-wave system, see || and the references 
therein. The crucial difference of (ITS) to the usual presentation of the 3- 
wave system is that we retain the linear terms. These terms usually are 
removed by the ansatz £ + ip^ = Aexp(it), etc., which leads to equations in 
the amplitudes A etc., see ||. In our treatment we retain the linear terms, 
because they determine the swing plane to lowest order. Keeping these terms 
enables us to find the swing plane without invoking the "pattern evocation 
in shape space" hypothesis, as was done in 0. 
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4 Reduction to one degree of freedom 

In what follows we study the geometry of the energy momentum map 

EM : T*R 3 -> R 3 : E i-> (#(E), J^S), J 2 (S)) (23) 

of the Liouville integrable system (if, J 1 , J 2 , T*R 3 , Q). This energy momen- 
tum map is related to the energy momentum map 



EM : T 



*TTT)3 



:Z^(H(Z),L(Z),H (Z)) 



of the resonant swing spring system by EM. = AoEMoty. Here \I/ is the linear 
symplectic map (14]) and A is the invertible linear map ~ 



2 \ ™ 

o -l i . lhus 

11/ 



the Liouville integrable systems (H, L, H , T*R 3 , u) and (H, J 1 , J 2 , T*R 3 , Q) 
are equivalent. 

The simplest way to reduce (H J J 1 J J 2 J T*R 3 ,u)) is to consider the two- 
torus action generated by the momenta J 1 and J 2 . Concretely, define a 
T 2 -action 



$ : T 2 x T*R 3 -> T*M 3 : ((t, s), 5) i-> <^ J1 o<p 



f(S) 
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(24) 



which comes from combining the S 1 actions (21) and (0). From their mo- 
mentum maps we obtain the momentum mapping 



J : T 



*TD>3 



>2 . = 



(J\E),J 2 (E)) 



(25) 



of $. To reduce the integrable system (if, J 1 , J 2 ,T*R 3 ,cD) by the symmetry 
$, we use invariant theory. In this approach, the generators of the algebra 
of invariant polynomials are used as new coordinates. The algebra of T 2 - 
invariant polynomials is generated by 



2 , /-2 2 i 2 2 i z-2 

Pi = P<: + £ , P2 = P v + V , P6=P( + ( 



(26) 



and 



p 4 = {£v-p&ri)pc-(&r,+'nPt)(> P5 = {€v-PtPv)(+(&v+vPt)Pc ■ (27) 
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The invariance of p± and ps is obvious if we write them as 
p 5 + ipA = (£ - ipz)(v ~ W V ){C + iPc), 

Note that the cubic part of H is exactly A p&. The invariants are subject to 
the relation 

pl + p\ = P1P2P6, pi > 0,p2 > 0,pe > 0. (28) 

Therefore the space Pj lt j 2 = J~ l {jii 32} /T 2 of orbits of $ with momentum 
(J13J2) is defined by fl2"8|) together with 

pi + p 6 = 2ji, p 2 + p 6 = 2j 2 , (29) 

which just expresses J 1 and J 2 in terms of invariants. Eliminating pi and p 2 
gives 

G(pa, p 5 , Pb) = P\ + P\ - Pe (2.7*1 - p 6 )(2j 2 - p 6 ) = , (30) 

where < p§ < min(2ji, 2j 2 ). This defines the fully reduced space Pji j2 as a 
submanifold of 1R 3 with coordinates (p4, ps, pe). Because of the restriction on 
p 6 and since pl + p\ is nonnegative, there are essentially three possibilities 
for the geometry of Pj lt j 2 , see Figure [[]. These possibilities are determined 
by the position of the roots of the polynomial 

P 3 (p 6 ) = p 6 (2j 1 -p 6 )(2j 2 -p 6 ). (31) 

1- ji 7^ 32 7^ 0. The space J _1 (ji, 32) is smooth and the action $ is free. 
Then the fully reduced space is either 

P h,h -Pl + Pl = Pe(2ji - p 6 )(2j 2 -pe), 31 > 32 > 0, < p 6 < 2j 2 , 

or 

P ^2 -Pl + Pl = Pe(2j'i - p 6 )(2j 2 - Pe), h > 3i > 0, < p 6 < 2j 1; 

which are each diffeomorphic to a smooth 2-sphere. Every point on the fully 
reduced space reconstructs to a two-torus orbit of the free action $ in phase 
space. 

2. jt = 32 7^ 0. The space J~ x {j\, ji) is smooth and the action $ has a fixed 
point (p 4 , p 5 , pe) = (0, 0, 2ji) with isotropy group S 1 . Then the fully reduced 
space 

P h,h -pI + pI = Pe(2ji - p 6 ) 2 , Ji > 0, < p 6 < 2j 1; 
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Figure 1: The polynomial P3 (upper row) and the corresponding reduced 
spaces Pj lt j 2 (lower row) illustrating the three types of reduced space using 
ji = 1, h = 3/2 (left), ji = j 2 = 1 (middle), and j x = 0, j 2 = 1 (right). 



is a topological 2-sphere with one (conical) singular point (0,0,2^). This 
singular point reconstructs to a pure springing motion on £ = rj = p^ = p v = 
and ( 2 + p"l = 2ji. The general point of the fully reduced space reconstructs 
to a two-torus in phase space which has angular momentum zero. This two- 
torus dynamically decomposes into an S l family of planar S 1 motions. 

3. ji = or j 2 = 0. The fully reduced spaces -Pj^o, ji > and Poj 2 , 
j 2 > each are a point = p$ = p§ = 0. In T*M 3 this reconstructs to a pure 
swinging motion with nonzero angular momentum on the circle £ 2 +£>^ = 
r] = p v = C — P( — (counterclockwise in the (x, y) plane projection) or the 
circle if + p^ = 2j 2 , (,=Pt = (=P( = (clockwise), respectively. If 
ji — h = 0, then P ,o is the point p 2 = P5 = Pe = 0, which reconstructs to 
the equilibrium point £ = p^ = 77 = p v = ( = p^ = 0. 

Using physical language, one would say that the projection of a generic 
motion of the resonant swing spring to the (x, y) plane is elliptically, linearly, 
or circularly polarized, corresponding to cases 1, 2, 3, respectively. 

Because the Hamiltonian H (17) is invariant under the T 2 action $, it 
induces the fully reduced Hamiltonian 

H juj2 : P jltj2 Ctf^K: (p 4 , p 5 , p 6 ) h-> d \ + j 2 + Ap 4 . (32) 
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Figure 2: The intersection of the singular reduced space P^i with the plane 
h = const gives an integral curves of the reduced one degree of freedom 
system. 

The integral curves of the reduced system are now given by the intersec- 
tion of the reduced phase spaces {G = 0} and the planes {Hj 1 j 2 = h}, as 
illustrated in Figure Q The fact that the Hamiltonian is a linear function 
does not imply that the system is a linear dynamical system, because the 
nonlinearity is contained in the Poisson bracket that gives the dynamics on 
the reduced space. A direct calculation shows that on C 00 ^ 3 ), where M 3 has 
coordinates (p4, p$, p$), the original canonical Poisson bracket on phase space 
(TM 3 , to) induces a Poisson bracket { , } on M 6 with coordinates (p\, . . . , p?,). 
Its nonzero brackets are given by 



These brackets have the generators of the symmetry pi + p 6 and P2 + p& 
as Casimirs, which we then use to eliminate pi and p2- This shows that a 
Poisson bracket on the fully reduced space Pj lt j 2 (considered as a subset of 



{P4,P5} 



{p2,pi} = {p4,Pe} = ~2p 5 , 
{p2,pn} = {p5,pe} = 2p 4 , 

(Pi + P2)PS - PlPl ■ 



(33) 
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IR 3 with coordinates (p 4 ,p 5 ,p 6 )) is 



dG 

{p6,pi} = W~ = 2p 5 

dG 

{P4, Ps} = tt- = Pe(2ji + 2j 2 - p 6 ) - (2j! - p 6 ) (2j 2 - p 6 ) (34) 
op 6 

{P5,Pe} = t^— = 2p 4 , 

OpA 



where G is given in (|30|). 

Hence the integral curves of the fully reduced Hamiltonian vector field 
X H jld2 on ( p iij 2 >{> }) satisfy 

p 4 = {P4,%jj = 

p 5 ={ps,H jl>j2 } =\P' 3 (p 6 ), see© (35) 
p 6 = {p 6 ,H jlj2 } = 2Ap 5 

When ji = ]2 > 0, the 2ji-level set of the fully reduced Hamiltonian Hj 1 j 1 is 
the intersection of the 2-plane {p 4 = 0} with the fully reduced space Pj 1 j 1 - 
This level set is an orbit of Xu- , which is homoclinic to the conical singular 
point (0,0, 2ji) of Pj 1 j 1 . The conical singular point reconstructs to the hy- 
perbolic periodic orbit 7 of pure springing motion; while the homoclinic loop 
reconstructs to a 2-torus bundle over each point of Hj^iT.]^ \ {(0, 0, 2^)}. 
Thus f/~j i (2ji) reconstructs to the stable and unstable manifold in T*R 3 of 
the hyperbolic periodic orbit 7. This invariant manifold is topologically the 
product of a once pinched 2-torus and a circle. 

Describing the geometry of the 3-tori and the dynamics of the swing 
spring near this homoclinic invariant manifold is the main objective of the 
remainder of this paper. 



5 The critical values of EAi 

We now determine the set of critical values of the energy momentum mapping 
SM (H) of the Liouville integrable system (#, J 1 , J 2 ,T*R 3 ,D). The set of 
critical values of EM. is called bifurcation diagram S. The critical values of 

EM. can be are determined from the equilibrium points of the reduced system. 
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Figure 3: Set of critical values of the energy momentum map (bifurcation 
diagram). It consists of two smooth patches intersecting in two lines and a 
thread (thick black line) connected to the former only at the origin. The 
image of the momentum map contains the thread, hence it is not simply 
connected. 

Geometrically this occurs for energies h for which the plane {Hj 1 j 2 = h} 
intersects the fully reduced space Pj lt j 2 in a (singular) point. This occurs 
in different ways for the three types of reduced spaces. In case 1, it occurs 
through a tangency of {Hj lt j 2 = h} with Pj lt j 2 - In case 2, it is either again 
a tangency or the plane {Hj 1 j 2 = h} contains the conical singular point 
(0, 0, 2jx). Hence {ji,j 2 , h) = (ji, ji, 2ji). In case 3, either j x or j 2 is zero, so 
that the critical values are (ji,0, ji) or (0, j 2 , j 2 ), respectively. 

The points of tangency between {Hj lt j 2 = h} and Pj lt j 2 flSTf ) can be ob- 
tained using Lagrange multipliers. This leads to the condition that the poly- 
nomial 

Q(pe) = ~(ji+h ~ hf + p 6 (2ji - p 6 )(2j 2 - pe) (36) 

has a multiple zero in [0, min(2j 1 , 2j 2 )]. When ji = j 2 > and h = 2j 1; we 
see that 2ji is a multiple root of Q in [0, 2ji\. 

The following argument determines a parametrization of the bifurcation 
diagram E. The bifurcation diagram S is contained in the set of (h,jx,j 2 ) G 
M 3 where the polynomial Q (|36D has a multiple complex root. This set is 
the discriminant locus of the polynomial Q. If we restrict ourselves to the 
set of {h,ji,j 2 ) G M 3 where the polynomial Q ( j36[) has a multiple root in 
[0, min(2ji, 2j 2 )], then we obtain the restricted discriminant locus A of Q, 
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which is equal the set of critical values E of the energy momentum mapping 
EM.. When (h, j\, 32) G A, then Q may be written as 

Q(Pe) = pl- 2(j'i + ]2)pl + 4j'ij 2 Pe - ^ 0*i + - /^) 2 

= (p 6 - s) 2 (p 6 - t), s G [0, min(2j 1 , 2j 2 )] 
= pl - (t + 2s)p 2 + (2st + s 2 )p 6 - s 2 t. 

Comparing coefficients gives 

t + 2s = 2(j 1 +j 2 ) 

2st + s 2 =4jU2 (37) 

s 2 t =&(ji+h-h)*. 



Since s > 0, from the third equation in (j37|) it follows that t > 0. Taking the 
square root of both sides of this equation gives 

3\ + 32 - h = e^sy/t, (38) 
where E\ — ±. Using the first two equations in (3"7p, we see that 

4(ji-j 2 ) 2 = 40 1 +j 2 ) 2 -16jU2 

= (t + 2s) 2 -4(2st + s 2 ) = t(t-4s), (39) 

which implies that 

either t > 4s and t > 0, or t = 0. (40) 

When either of the conditions ( f40l) hold, we may take the square root of both 
sides of (139|). We then obtain 



ji -h = \e 2 ^Jt{t-4s), (41) 
where £2 = ±. Adding and subtracting (j38"D from 

3i+32 = \(t + 2s) } (42) 
(see the first equation in ([37])) gives 



j! =3i(s,t) = \(t + 2s + £ 2 y/t(t-As)) 

32 =h{s,t) = \(t + 2s-£ 2 y/t{t-4 S )). 



(43) 
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Subtracting ( f4"2|) from ( pl£| ) we obtain 

h = h(s, t) = \ (t + 2s) - £ X A sv^. (44) 

Consider the mapping 

V £l , £2 : D C R 2 -> M 3 : (s,t) h-> (A( s ,t), j 2 M)), 

where -D is the closed subset of IR 2 which is the union of the intersection of 
the closed half planes {s > 0}, {t > 0}, {£ — 4s > 0} and the closed half line 
L_ = {(s, 0) G M 2 | s > 0}. For every choice of sign of e\ and E2, the map 
V £l;£2 parametrizes a patch of the bifurcation diagram E. All the patches 
together are shown in Figure [| 

The image under V £l)£2 of the closed half line L_ is the thread X = 
{ (2ji, jijji) G IR. 3 1 ji > 0} , which is attached to the 2-dimensional pieces 
V £l , £2 (D\L_) of E only at the origin (0, 0, 0). Since the range X of the energy 
momentum mapping SAi (^3]) is the closed region in M 3 bounded by E and 
containing the thread T, we see that the set of regular values X \ T in the 
image of EM. is not simply connected. Thus the Liouville integrable system 
(H, J 1 , J 2 , T*M. 3 , lo) possibly can have monodromy. The thread T represents 
the unstable springing motion including its separatrix. The intersections of 
two patches, which form 2-dimensional pieces in the boundary of the image 
X, represent the periodic orbits of left or right circular swinging motion. Any 
other point on the boundary represents two-tori. Their special feature is that 
the instantaneous ellipse in the (x, y)-plane formed by their projection has 
constant excentricity, see section 8 below. All the remaining points in the 
image of £ Ai represent generic motion on three-tori. 

For the geometrically inclined reader, we will show in sections 6 and 7 
that the integrable system (H, J 1 , J 2 , TM 3 , cD) has monodromy. Analysts may 
skip these sections on the first reading and go directly to the derivation of 
the swing plane angle in section 8. This provides an analytic proof of the 
fact that the resonant swing spring has monodromy. 



6 A different reduction 

To uncover the geometry of the energy momentum map EM. (J23|), we reduce 
the Liouville integrable system (H, J 1 , J 2 , T*M 3 , u) to a two degree of freedom 
Hamiltonian system by removing the ^-symmetry (pj (21). As removing 
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the symmetry induced by (p s fl2"2"|) leads to the one degree of freedom system 
we have already discussed in section 4, we will omit this second reduction. 

To reduce the ^-symmetry ipf we use invariant theory. The algebra of 
polynomials on T*IR 3 , which is invariant under this symmetry, is generated 
by 

°\ = p\ + £ 2 °"3 = v o 5 = £pc - Cvi 

(45) 

0"2 = p} + C 2 4 = Pr, 6 = + PiPC- 

These invariants are subject to the relation 

a\ + a\ = a x a 2 , a x > 0, o 2 > 0, (46) 

which define the space T*M?/S l of orbits of the S^-action ip{ . The reduced 
space Pj 1 = ( J^^di) / S 1 of orbits of ip J on the jidevel set of the momentum 
J 1 (^) is defined by (|46| ) together with 

I ( ai + a 2 ) = 3l > 0. (47) 

Eliminating <j\ from fl46|) using ([|7]) gives 

a 2 + a 2 = o 2 (2jx - <j 2 ) , < a 2 < 2 , (48) 

which defines Pj l as a submanifold of 1R 5 (with coordinates (o 2 , . . . , cx 6 )). 
Since the Hamiltonian H (|T7D is invariant under yj/ 1 , it induces the Hamil- 
tonian 

Hji '■ P h ((T 2 ,(T 3 ,... ,<T 6 ) ^•il + |(0" 2 + Cr 2 +(T|)+A((73(T 5 -(T4(T 6 ), (49) 

where A = 

Consider IR 6 with coordinates (ax, o 2 , . . . , ae). On C°°(1R 6 ) there is a Pois- 
son bracket { , }, induced from the Poisson bracket on the space of smooth 
functions on (T*IR 3 ,cD), such that the nonzero brackets are 

{05,0-1} = 2o 6 {o 6 ,Oi} = -2a 5 {05,02} = -2o 6 

(50) 

{06,02} = 2o 5 {04,0-3} =-1 {o 6 ,o 5 } =oi-o 2 . 

Note that 01 + o 2 and o| + o| — oio 2 are Casimirs for the Poisson algebra 
(C°°(IR 6 ), { , },•)• The reduced Hamiltonian vector field X s . on R 6 has 
integral curves which satisfy 

01 = -2A(o 3 o 6 + o 4 o 5 ) d 4 =-o 3 -Ao 5 

2 = 2A(o 3 o 6 + o 4 o 5 ) 05 = -o 6 + Ao 4 (oi - o 2 ) (51) 

3 = o 4 — Ao 6 d 6 = o 5 + Ao 3 (oi — o 2 ). 
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7 Monodromy 



To show that the Liouville integrable system (if, J 1 , J 2 , T*M 3 , Q) has mon- 
odromy, we fix ji > and reduce the S 1 action tpj (^1|). We obtain the 
two degree of freedom Liouville integrable system (H^, J 2 1 Pj 11 { , }), where 
J 2 = \ (02 + + cr|) is induced from J 2 (^). We will show that this system 
satisfies the hypotheses of the monodromy theorem as stated in Cushman 
and Duistermaat, (see also Matveev M and Zung [|12| ). Thus it has mon- 
odromy. This shows that the original Liouville integrable system also has 
monodromy. 

First, we verify that the reduced space Pj t , j\ > is a smooth manifold. 
Since Pj 1 C M 5 is defined by 

F(a) = a\ + o\ - 2 3l a 2 + o\ = 0, < a 2 < 2j u 

the point a = (<r 2 , . . . , cr 6 ) is singular if 

(0,0, 0,0,0) = DF{a) = (2{a 2 - ji), 0, 0, 2<r 5 , 2a 6 ), 

that is, <t 2 = ji and cr 5 = <r 6 = 0. But F(j 1 , a 3 , cr 4 , 0, 0) = — jf ^ 0. Therefore 
Pj 1 has no singular points. 

Second, observe that the S 1 action on Pj x induced by the S 1 action <pj 
flU) is 
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— sin s 
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a. 



The action ^ has a unique fixed point cr° = (2j 1; 0, 0, 0, 0) on Pj 1 C M 5 . 
Since the reduced Hamiltonian 

H h = il + | (°2 + °3 + °4) + -M°"30"5 - <T 4 flr 6 ) 

is invariant under ^ , the point a is a critical point of Hj 1 . Because 
{Hj 1 , J 2 } = 0, the Hamiltonian vector fields Xq and Xj 2 commute. From 
([51]) and the fact that <Ji + a 2 = 2ji , we know that Xq on 



p5 ; 



IS 



2X(a 3 a e + a^a^)-^ + Xa, 



+ (a 6 - 2Aa 4 (ii - o-2))g|- - (a s + 2Xa 3 (ji - a 2 ))^-. 
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Clearly (a ) = 0. Moreover, the linearization of on the tangent 
space 



1 n 



TsoP h = ker DF(a°) = ker(2ji, 0, 0, 0, 0) = span { 



a a a a 



IS 



^ = DXft (a )\T^oP jl 



/ A \ 

AO 

2Aji 1 

V 2Aji -10/ 



Since the characteristic polynomial of Y is (x — Ay / 2j^") 2 (x + Ay / 2jj~) 2 , the 
equilibrium point a is hyperbolic. 

We now look at the energy momentum mapping 



£M n : P n C 



d2 . 



0'i + i (o"2 + o-f + + A(a 3 o- 5 - a A a 6 ), \ (cr 2 + of + af )) 



(52) 



of the Liouville integrable system (Hj 1} J 2 , Pj 1} { , }) with fixed ji > 0. The 
fiber SAij^h, j 2 ) is compact. To see this note that cr 2 > on Pj 1 and 
cr 2 + o"| + <t| = j 2 > 0. Therefore < <j 2 < j 2 , < a 3 < -y/72 an d < cr 4 < 
Vj2- But 



OV 



+ a 6 = 0-2(27'! - o- 2 ), < o- 2 < 2ji 



is the defining equation of P 3 - r Thus o" 5 and a 6 are bounded. Hence £Ai~^(h, 
J2) is compact. In other words, ^.M^ is a proper map. To show that 
SAij^h, J2) is connected, note that it is the total space of an S 1 bundle 
over the h- level set of the fully reduced Hamiltonian Hj lt j 2 ( p2|) on the fully 
reduced space Pj 1 j 2 (PB]). Since this level set is connected, when (h,ji,j 2 ) lies 
in the image of EM. (p3|) (and thus (h,^) lies in the image of SAi^ ([52])), 
it follows that £AiJ 1 1 (h,j 2 ) is connected. At the critical value of 
£M.j 1 we have already observed that, reconstructing the critical 2ji-level set 
of the fully reduced Hamiltonian Hj 1 j 1 on Pj 1 j 1 , we obtain a once pinched 
2-torus, which is the critical fiber EAiJ^ftji, ji). Thus we may apply the 
monodromy theorem. We find that over a loop T in the set of regular val- 
ues in the image of 8Atj 1 , where T encircles the critical value (2ji,ji), the 
2-torus bundle 8 Ai ^ (r) has monodromy (\ M . 
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8 The swing plane angle 

In the original phase space T*R 3 with canonical coordinates (x, y, z, p x , p y , p z ) 
let us project the motion of the swing spring onto the (x, y) plane of its 
configuration space M 3 . Taking into account only the relevant quadratic 
terms 

H° xy = l(pl+p 2 y + x 2 + y 2 ), (53) 

of the Hamiltonian H ( |TTD of the resonant swing spring, we obtain the Hamil- 
tonian of a two degree of freedom harmonic oscillator. The projected motion 
of this oscillator is an ellipse E. For the interpretation of the solutions of 
the resonant swing spring we need some elementary facts about these ellipses 
which are the content of the following 

Lemma 1. The center of E is at the origin and its major axis makes an 
angle 

# = I tan' 1 2 ( x y+P*Pv) (54) 
2 p x + x 2 — p 2 — y 2 

with the x-axis, assuming that p 2 + x 2 > p 2 + y 2 . Otherwise, d is the angle 
between the major axis and the y-axis. In the coordinates (£, 77, CjP£jPj?>Pc) 
on (T*R 3 ,cD) this angle is 



# = W = Itan- 1 ;' " V = --itan-^ + itan- 1 ^. (55) 
tPv-VPi 2 Pa Pv 

The area of the ellipse E is given by 7r(xp y — yp x ) 2 = ttL 2 and its eccentricity 



is 



e(E) = 2(K + K- 1 y 1 , k 4 = ^^. (56) 



Proof. Recall that (see [§) 



n = x 2 + p 2 x r 3 =xy + p x p y 

t 2 =y 2 + pi n = xp y - yp 



(57) 



■y 

are integrals of the 2-dimensional harmonic oscillator vector field 



X - ® + ® x ® ® (58) 

H xy P x Q x PV y X y ■ 

Moreover they satisfy the relation 

rl + tI = nr 2 , n > 0, r 2 > 0. (59) 
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The projection of a motion of the harmonic oscillator with initial condition 
(x°, y ,p®,Py) and positive energy onto the (x,y)-plane is the curve 

7 : 1 1— ► 7(2) = (x(t), y(t)) = (x cot t + p x sin t, y cost + p y sin t). (60) 

The initial condition determines the values of the integrals Tj, say rf. Then 

- x(t) 2 )(r 2 ° - y(t) 2 ) = ( Px (t)Py(t)) 2 = (r 3 ° - x{t)y{t))\ 

which after some simplification gives 

r° 2 x(t) 2 - 2r° 3 x(t)y(t) + r^t) 2 = - (r 3 ) 2 = (r 4 ) 2 . 

Since t° + r° = > and t^t® — (rf ) 2 > 0, the quadratic form 



(%,y) 



0\2 



(61) 



is positive definite. Hence the image of the curve 7 is an ellipse E. 

Changing coordinates by ( M = ^ ™ ( 8 f M transforms the quadratic 

form Q (^n|) into the quadratic form 



C 2 TS + S 2 T° 



l -2cstI 



cs(t° 2 - t 2 °) - (c 2 - s 2 )r ; 



2 wo 
3 



C S (Tf-T 2 )-(c 2 - S 2 )T 3 

2csro° 



2 ^0 1 „2 ^0 



S TK + C T- 



1 "T ^<-f '3 
^ ^0 



(62) 



where c = cos 9 and s = sin#. Suppose that rf > tj. Choose 6* so that 

2r° 

= Uan' 1 3 (63) 



Then the quadratic form fl5"2] ) becomes A_ £ 2 + X + rf = (r^) 2 , where A± are 
eigenvalues of Q, that is, 



A 



± 



(r°! + r° 2 ) ± yj (t? + t§)* - AtfY* 



> 0. 



(64) 



Since A + > A_ > 0, a piece of the £-axis is the major axis of E and 9 ( |63| ) is 
the angle between the major axis and the x-axis. Note that the eccentricity 
e of the ellipse E is e 2 = 1 



A_ 
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Figure 4: The eccentricity e(pe) for ji = 1, j'2 = 1-5. The graph is even closer 
to one when ji and j'2 are closer to each other. 



We now compute k. By definition 

2ji - p 6 Pi £ a + P\ 



So 



K 2 



Therefore 



2j 2 - Pe P2 ^ 2 + P 2 

x 2 +p 2 x + y 2 +p 2 y - 2(xp y - yp x ) T j_ + r 2 ° - 2r 4 ° 
+ p 2 +y 2 +p 2 + 2(xp tf - j/p x ) ~ r 1 ° + r 2 + 2r 4 ° 



x 



(rf + r°) 2 - 2(r°) 2 
( t o + t o + 2t O)2 • 

Tf + r 2 ° + 2r 4 ° 

, since the characteristic polynomial 

A+ + A_ + 2v/A+A_ of Q ig A 2 _ (r o + r o )A + (r o )2 

1 - 



2 _K (1 - k 2 ) 2 

A_ (1 + K 2 ) 2 ' 



e 



which implies that e 2 = T , ^, k m , that is, e = = , 2 _i ■ □ 

For small but non-zero angular momentum I the motion goes from a small 
circle to a very eccentric ellipse, see Figure [|. When I = then j\ = H, hence 
K — 1 and e = 1: the motion takes plane on a line, always. However, when 
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ji 7^ 32 then the eccentricity will move between and 1 — 21/ j. E is close to 1 
for most values of p 6 , and $ gives the angle of the swing plane. Only when p 6 
is near its maximum the eccentricity becomes small. For the full Hamiltonian 
H of the swing spring, the angle •& is no longer constant, but instead is slowly 
varying. In this situation, *& is the instantaneous angle of the swing plane. 
Similarly p 6 is changing in time, so e becomes the instantaneous eccentricity. 
Even in the full dynamics the instantaneous ellipse E has constant area 7rL 2 , 
because L is the conserved angular momentum. In general the eccentricity of 
E changes during the motion of the swing spring. However, since e(E) only 
depends on p 6 , we see that it is constant for the critical 2-tori corresponding 
to a stable equilibrium point of the fully reduced vector field . . The 
swing angle $ should not be confused with the polar angle of polar coordinates 
in the (x, y)-plane, even though •d is canonically conjugate to L. As we will 
see the angle i? commutes with H , which is not true for the polar angle. 

We now construct the action-angle variables corresponding to the T 2 sym- 
metry group. In section 3, we have found a Hamiltonian 2-torus action ([24]) 
on T*M 3 with momentum map J = (J 1 , J 2 ) (|25|). For a regular value {31,32) 
of J, the two-torus action on the level set 32) = (^ 1 )~ 1 (ji)n(J 2 )~ 1 (j2) 

is free. Hence we obtain a smooth 2-torus bundle 7rj lj2 : J _1 (ji, j 2 ) — > Pj u j 2 
over the fully reduced space Pj lt j 2 . We would like to find action- angle co- 
ordinates for the 2-torus fiber F p = irj x ,- 3 (p) over the point p E Pj 1 ,j 2 - The 
Hamiltonian vector fields Xji and Xj2 on (T*M 3 ,cD) leave each 2-torus F p 
invariant and have only periodic orbits of period 2n on F p . The same is 
therefore true for Hq = J 1 + J 2 and L = J 2 — J 1 . Since we are interested in 
the angle 9, we use H and L as action variables instead of J 1 and J 2 . We 
denote the values of these actions by j and I, respectively. This leads to the 
following 



Lemma 2. For all regular values of J the variables (ii, 9\, I2, ^2) given 
by h = H , I 2 = L, (see 4M) and {jffl)), 



9 1 = ^tan~ 1 —, 9 2 = \ tan" 1 — - \ tan" 1 — (65) 
PC P v Vi 

and p4,P5,pe on Pj 1 j 2 are a coordinate system for an open subset of phase 
space T*M 3 . In these coordinates, the Hamiltonian of the resonant swing 
spring is 

H = h + \p 4 . (66) 
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Moreover, the following Poisson bracket relations hold. 

(67a) 
(67b) 

(67c) 



{0i,h} = 


1, 


{^2,/ 2 } = 


1, 




1 P4 

2 ' 
P6 


{^P 5 } = 


1 P5 

2 ' 
P6 


{#2,P 4 } = 


P 4 -^2 


{#2,P 5 } = 


P5^2 
P1P2' 


{Ol, P6 } = 


1, 


{p 4 ,P 5 } = 


Pe(Pi +P2) -P1P2, 


{p5,pe} = 


2p 4 , 


{p 4 ,Pe} = 


-2p 5 , 



(67d) 
(67e) 

where p\ — I\ — l 2 — Pe and p 2 — h + I 2 — p§- All other brackets vanish. 

Proof. This last assertion can be verified by direct calculation using the old 
variables (£, 77, (, p v , p^) . For example 

{02, h] = (pc^-C4)tan-^ = I. 
h} = 

{02,h} = (+^+^-£4-^)^2 = 

{9 2 ,i 2 } = (- p ^ +Pn f q+ ^-^e 2 = i 

- Pt;Pv)P( - (&v + m>z)t _ P4 



1 _ 1 

2 2 



{0!, P4> = 
{^2,P4} = 



2(p 2 + C 2 ) 2p 6 
(p| + £ 2 - P% - V 2 )((£V ~ P&v)Pc - (&v + VPt)0 __ PaXt. 



2(p 2 +e)(P 2 +^ 2 ) Pip 2 ; 



and similarly for the other brackets. That all brackets between the new 
variables can be expressed in terms of the pi follows from the fact that 

1) 1 1 and I 2 are actions; 

2) the flow of the Hamiltonian vector fields of I x and I 2 leaves pi invariant; 

3) the actions I\ and I 2 and the angles 6±, 9 2 form a symplectic coordinate 
system in the two-torus fibre of the bundle 7ij lt j 2 . □ 

Lemma 2 allows us to give a description of the motion of the resonant 
swing spring in phase space in such a way that the fully reduced system is 
an invariant subsystem that is driving the dynamics in the fibres. Using the 



22 



fact that Ap4 = h — j, the equations of motion for the angles in the fibre are 

e 1 = {9 1 ,H} = l + f^-, (68a) 

02 = {e 2 , H} = {h - 3)l . (68b) 

[P6 - 3 ~ 1){P6 ~ 3 + I) 

They are driven by the solution of the the second order differential equation 

p 6 = 2A 2 P^(p 6 ) = 2A 2 (2p 6 (j-p 6 )-(p 6 -j-/)(p 6 -j + /)) (69) 

Equations ( |68b| ) and fl6"9|) can be found in |J. We rederived them because 
in [||] they were derived under the hypothesis of "pattern evocation in shape 
space", which might have been an approximation. From our derivation we 
see that no approximation is involved and that the angle of the swing plane 
is simply an angle of an action-angle coordinate system. 



9 Actions and Rotation Numbers 

We are interested in the change of # 2 over one period of the motion of p§. 
First we directly derive an integral for this change A# 2 by changing the time 
parametrization in ( |68b|) . Then we will show that A# 2 is a rotation number. 

The integral curve in the fully reduced system is given by the /i-level set 
of the fully reduced Hamiltonian Hj lt j 2 ( |66"D on the fully reduced space Pj u j 2 
(p0|). Hence the motion takes place on a family of real affine elliptic curves 
£h,ji,j 2 defined by 

p 2 = Q(p 6 ) = -— (h - j) 2 + p 6 (p 6 - j - 0(pb -3 + 1), (70) 

when < p6 < min(2ji, 2j 2 ). Equation (|70|) is obtained by eliminating p4 
from ( P3|) using h = ji + j 2 + Ap4 (which defines Hj^ 2 {h)). This exactly 
reproduces the polynomial Q already defined in (|36|). Therefore also the 
analysis of the double roots done for the energy momentum map also applies 
here. The equation (|69| ) for p 6 can be integrated once and we find 

p 6 = 2Ap 5 = 2A V / Q(p 6 ) • (71) 

Note that this is one of the fully reduced equations of motion (|35|). When 
(^ ji) J2) is a regular value of the energy momentum mapping EM. (|23|) of 
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the swing spring (which we henceforth assume), then the £h,j u j 2 * s smooth. 
When h = j the polynomial P$(pq) = Pe(j ' + 1 — pe) (j ' — I — Pe) is nonnegative 
on [0, min(2ji, 2j 2 )]. Therefore the polynomial Q has three distinct real roots 

< p 6 < pt < min(2j 1 , 2j 2 ) < p° 6 . 

The motion of the fully reduced system on if _1 (/i) takes place when p 6 lies in 
[pg ,Pg]. Since Hj^- 2 {h) is diffeomorphic to a circle, when the fully reduced 
motion runs through a period the time parameter p§ traverses the interval 
[pg , pg] forward and backward once. Accordingly the period of the driving 
motion is given by 

T = 2 j"" i^- (72) 
J P - 2Ap 5 1 ; 



To find the change of 8 2 over this period, we introduce in (|68a|) and ([38 
a new time scale p6 defined by fl7j} ). We obtain 

dfc = (l + ^fA 4^ (73a) 

\ 2 ft j 2Ap 5 1 ; 

= 7 i {k Z j)l , ^ (73b) 

which are differential forms on £h,ji,j 2 - Therefore the change in the swing 
angle during a period of the motion of the swing spring is 



A6 2 (h,j,l) = 2 f 6 dd 2 

(h-j)l 1 d p 6 



A J p - (j + l- Pe) (j -I- Pe) ^Q{pa)' 



(74) 



To see that A8 2 is in fact a rotation number of the integrable system 
(H, J 1 , J 2 ,T*M. 3 ,Q), we have to compute the third action. First we find the 
angle # 3 conjugate to p§ on Pj 1: j 2 \ {(pg ,0)}. A calculation shows that 



^ 3 = |tan- 1 ^. (75) 

P4 

does the job, because {83, p$} = 1. Therefore, we can define the third action 
J3 as the integral of the canonical one form p 6 d8 3 over ^j b j 2 , namely, 

= ^ u** = *=i r + (3- 2A - &\ (76) 



2tt / 8ttA J p - V Pi P2 / P5 



24 



The last equality above is verified as follows. Computing the derivative of 
#3 along the integral curves of the fully reduced vector field Xj u j 2 and using 
(HI) gives 

d6> 3 x p 4 dp B h-j A 



dt 2 pl + pi dt 2A P1P2P6 

So 

d0 3 h-j {p x p 2 - p 2 p e - pipe 



PsiPe). (77) 
(78) 



dp 6 4A P1P2P5P6 
Substituting p 6 = 2j 1 — pi and p 6 = 2j 2 — p 2 gives the equality (|76|) . 

The integral ([TJJ) implicitly defines the Hamiltonian if as a function of the 
three actions h, I 2 , I3. Hence we can obtain the frequencies of the conjugate 
angles 61, 6 2 , and #3 by implicit differentiation. In particular we have 

Rewritten this gives 

dH = -{-dh) w dIl -{-dh) -dr dh + {-dh) dh 

dH Jr <9# Jr <9# Jr 
= 77— dh + — - dl 2 + — - d/ 3 . 
dh dl 2 dh 



(80) 



(81) 



Hence 

■ = dH = fdh 
3 dl 3 V dh 

The above expression yields the period 2n which is the same as that given 
in ([721). To see that the integral (|72| ) and the one obtained from ( |8T| ) give 
the same result, one has to add a total differential to the integrand of the 
loop integral ([T2|), which does not contribute to the period. More precisely 
the identity 

!!* » (*zi( 3 -^-^)**)+d(*), (82) 
p 5 9/1 V. 2 V Pi P2 J P5 J VPs/ 

holds on £ h j 1 j 2 , where p 5 is a function of p 6 given by (f70|) . The frequencies 
of the other angles are similarly obtained as 

. 0if _ dh/di 

02 ~ df 2 ~ ~dhldh (83) 
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Figure 5: Contour plot of x f° r j — 1 5 I — 0.15. The contours 0, —0.05, —0.1 
are shown. 



and similarly for 0\. Hence we obtain the rotation numbers 

Now we show that A8 2 = — W23. The equality holds because the differentials 
only differ by a total differential, which does not contribute under the closed 
loop integral. By direct calculation one can check that 

d92 = a/(^r( 3 -7r^j7rJ + T 6 W' (85) 

where d(9 2 is given by ( |73 b| ) . The first derivation of an integral for the stepwise 
precession of the swing plane is simpler than the present one. But the latter 
derivation shows that the stepwise precession of the swing plane is one of the 
two rotation numbers of the invariant tori of the integrable approximation. 

In order to understand the solution curves of the differential equation for 
d9 2 / dp 6 ( |73b|) we can plot its solution curves. This is easily done by finding 
a function of 8 2 and p$ that is constant on the solutions. Such a function can 
be defined as 

y 2 
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where x(p 6 ,d0 2 / dp 6 ) = —{h — j) 2 /X 2 . A contour plot of x is shown in 
Figure || On the right part of figure the p 6 is large, so that the eccentricity 
is small. On the left side p§ is small, hence the eccentricity is close to 1, see 
Figure |[ 

10 Analysis of the swing angle 

In this section we examine more closely our formula (|74|) for the swing angle 
of the swing spring. We prove the following 

Proposition. The rotation number describing the angle of stepwise preces- 
sion for the integrable approximation of the resonant swing spring is given 
by 

A^ 2 = tan' 1 - ^ + O (yp + ((h-j)/\)*/j/j) . (87) 

Equation ( jS7f) shows that the swing angle A9 2 is a multivalued function 
of the parameters h, ji and j 2 - From this one can read off that the Liou- 
ville integrable system (H,L,Ho,T*M. 3 ,uj) describing the swing spring has 
monodromy. 

Proof. In order to simplify our analytical study of (f74[) we remove an un- 
needed parameter by the rescaling 

p 6 =jz, I = jb and af /2 = (88) 
Equation (173) then becomes 

A9 2 = ab f ^ dZ (89) 

J 2- (1 - ~Z) 2 ~ b 2 _ ~y _ 6 2) _ S 2 

Here < z~ < z + < 1 < z° are distinct nonnegative roots of the polynomial 

z((l-z 2 )-b 2 )-a 2 . 

Next we restrict the parameters a and b in (^) to lie on a line through 
the origin in parameter space M? and assume that they are both small. In 
other words, we introduce a small parameter e such that 

a = ea and b = eb, (90) 
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where (a, b) are fixed in parameter space. 

Now we find the Taylor expansion of the roots z~ , z + , and 5°. A calcu- 
lation gives 

r = aV + (2a 4 + aV) £ 4 + 0( £ 6 ) 
z+ = 1- Va 2 + 6 2 e - ± a 2 e 2 + 0(£ 3 ) (91) 
2° 



1 + Va 2 + o 2 e + | a 2 e 2 + 0(e 3 ). 

Because the root z + and the pole at 1 — eb coalesc at 1 as e — > 0, we 
introduce a shifted, scaled, and inverted new variable z by 

z=l + eb/z. (92) 

The inversion ensures that the new integration boundaries are finite in the 
limit e — > 0. Then (189T) becomes 



A ^ 2 = / ; „ , . „ I ,, . ,„ , ttt ' (93) 



az 2 dz 
^ - 1 v/-^[(a 2 + 6 2 )^ 3 - 6 2 z + eb 3 (z 2 - 1)] 

Factor the polynomial under the square root as 

-z[(a 2 + b 2 )z 3 - b 2 z + eb\z 2 - 1)] = (z - z~){z + - z)z(z - z°)(a 2 + b 2 ) 

and introduce a new integration variable by 

2^ = (2+ + z - ) + (z+-2T) cos 0. (94) 

The integral Q93|) becomes 

A# 2 = ^=L== I* f 44, (95) 



Va 2 + & 2 io 

where 



az 2 



^ (2 2 - 1)^(2 -2°) ' (96) 



Using (|92|) we find that fl9~i~D becomes 



= , e — —e + U(e ) 

Va 2 + 6 2 a 2 + 6 2 v ; 

z+ = -be + 0(e 3 ) (97) 
z° = + / b e - + 0(£ 2 ). 
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(98) 



(4a 2 + (3 - cos <p) (1 + cos 0)6 2 ) ^(cos0 - l)(cos</> - 3) ' 

Note that the error term is uniformly bounded in the interval of integration. 
The main purpose of the above transformations was to achieve this. Now the 
zero order contribution in (|95|) can be calculated. The substitution cos0 = 
1+2 sin-0 removes the root and the second substitution cosip = x rationalizes 
(|98|). Hence the integral becomes 



Using the above method, one can actually compute one more order. But 
it turns out to be zero. At 0(e 2 ) the integrand is not uniformly bounded. So 
more sophisticated methods are needed to obtain the first nonzero correction. 
Note that the error term is small when we are close to the thread in the 
bifurcation diagram. Then I and h — j are close to zero. Of course j itself 
is also small, since we must be close to the equilibrium, but I and h — j are 
considered to be much smaller. 

The result means that for each fixed value of j > 0, the swing angle A6 2 
is a multivalued function of I and h. This is another proof of the fact that 
the system has monodromy. The resonant swing spring provides an example 
in which the monodromy can be easily observed. 
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(99) 



Undoing the scaling gives the desired result. 



□ 
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